Phase coherence times in the multiple scattering of photons by cold atoms 
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We present an analysis of the dephasing present in the 
multiple scattering of photons by atoms with a quantum in- 
ternal structure. The corresponding phase coherence times 
T0 are obtained as a function of the Zeeman degeneracy of 
the atomic dipole transition and the polarization state of the 
photons. These results allow for an explanation of the recent 
experiments on coherent backscattering of photons from a gas 
of cold rubidium atoms where the height of the backscattering 
cone depends on the atomic internal degrees of freedom cou- 
pled to the polarization of the photons. Some consequences 
of these results are presented, and analogies with the case of 
electronic systems are highlighted. 
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Multiple scattering of waves in complex media ac- 
counts for a large number of coherent effects which are 
often presented under the shortcut of mesoscopic physics. 
These coherent effects result from interferences between 
the amplitudes associated to multiple scattering paths of 
the wave. In the weak scattering limit for which kl e 3> 1, 
where k is the wave-vector and l e the elastic mean free 
path (the average distance between successive collisions), 
the relevant physical quantities can be obtained from the 
probability P(r, r', t) of quantum diffusion defined as the 
disorder averaged probability to propagate a wave packet 
between the points r and r' in a time t [[l]j|] . It allows to 
describe either weak localization effects in metals jl],[|, 
spectral quantities or the coherent albedo and dynamical 
effects in multiple scattering of light by classical scatter- 
ers|. 

Interference effects are very sensitive to dephasing. 
Roughly speaking, a dephasing may originate either from 
an external field or from additional degrees of free- 
dom which affect in different ways each of the interfering 
amplitudes. An example is provided by the spin-flip scat- 
tering in metals where the spin of the electron rotates due 
to scattering by magnetic impurities || . 

In the presence of dephasing, the probability of quan- 
tum diffusion can be written as P(r,r',t) (e 1 ^*- 1 ), where 
the random variable <f>(t) is the relative phase of the two 
interfering paths. Its distribution depends on the origin 
of the dephasing and we denote by (...) the average over 
this distribution. In most cases we have (e^W) ~ e"*/^ 
at least for long enough times t. The characteristic time 
T0 is the phase coherence or dephasing time. An expo- 
nential decrease of the probability of quantum diffusion 
does not necessarily describe a dephasing process. For 



instance, the intensity of an electromagnetic wave which 
propagates in an absorbing medium decreases exponen- 
tially. But this is not a dephasing process since it affects 
equally both the coherent and incoherent contributions 
by a decrease of the overall intensity. We propose to 
speak of dephasing only when the coherent part is re- 
duced with respect to the incoherent part. 

Recently, the dephasing induced by the Coulomb inter- 
actions in metals has been reconsidered in detail in view 
of the recent experiments on the dephasing time T^iT) 
of the electrons at the Fermi level where an unexpected 
saturation at low temperature has been observed || . 

The propagation of photons in complex systems ad- 
dresses similar questions and provides new sources of de- 
phasing like the motion of the scatterers || . It also pro- 
vides access to the type of scattering experienced by the 
photons (Rayleigh, Mie, resonant, Raman etc.) so that 
we have at our disposal both a theoretical description of 
the elementary scattering events and its generalization to 
multiple scattering. 

Dephasing and decoherence in multiple scattering are 
sometimes perceived as a nuisance which prevents the 
observation of full interference effects. Instead, when the 
dephasing is well controlled, it should better be viewed 
as a unique tool in order to obtain new insights on the 
system that are usually out of reach in the very dilute 
and single scattering limits. A beautiful implementation 
of this has been demonstrated in the so-called diffusive 
wave spectroscopy 0. There, the controlled dephasing 
between multiple scattering trajectories of the photons 
resulting from the dynamics of the scatterers allows to 
probe this dynamics on time scales unreachable other- 
wise. 

Recently, the problem of multiple scattering of photons 
in a gas of cold atoms has been investigated in detail and 
coherent backscattering has been observed [00] . It has 
been shown that the existence of internal atomic degrees 
of freedom modifies significantly the coherent backscat- 
tering 10 1. The purpose of this letter is to present an 
analysis of the dephasing induced by the internal atomic 
degrees of freedom and to derive corresponding expres- 
sions for the phase coherent times as a function of the 
polarization state of the photons and of the Zeeman de- 
generacy. We believe that the present analysis provides 
relevant spectroscopic tools in order to probe the dynam- 
ics of cold atoms. 

We describe a gas of N atoms as two-level systems of 
characteristic transition frequency uo JlQ] . The ground 
state defines the zero of energy and has total angular 
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momentum J. The excited state has a total angular mo- 
mentum J e and a natural width T due to coupling to the 
vacuum fluctuations. We shall assume, moreover, that 
the velocity v of the atoms is small compared to T/k (k 
is the light wave- vector) but large compared to hk/M 
(M being the mass of the atom), so that it is possible to 
neglect the Doppler and recoil effects. The external de- 
grees of freedom of the atoms are therefore the classical 
assigned positions r a (a = 1, . . . , N) uncorrelated with 
one another. The corresponding one-atom Hamiltonian 
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describes the internal quantum degrees of freedom, in 
units where h = c = 1. No magnetic field is supposed to 
be present, so that the two levels are respectively (2J+1) 
and (2J e + l)-fold degenerate. 

The atom-photon interaction is described within the 
dipole approximation, namely using the Hamiltonian 
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where D a is the atomic dipole operator and E(r) is the 
quantized electric field operator (quantization volume V) 
evaluated at the center of mass r a of each atom. The nat- 
ural width of the excited atomic state is T — d 2 ujQ/3ireo 
where d = (J e ||D|| J)/s/2Jl + 1 in standard notations. 
We define the dimensionless dipole operator d = D/cf 
with non vanishing matrix elements only between the 
two states J and J e . The elastic scattering process be- 
tween the two states |ke, Jm) and |k'e', Jm!), where |ke) 
is a one-photon Fock state of the free transverse elec- 
tromagnetic field in the mode k of polarization e, is 
described by the single scattering transition amplitude 
tij(m,m' ,u>) — t{uj){Jm'\didj\Jm) . Here, the resonant 
scattering amplitude is given by 
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where 5 = u) — u>o is the detuning of the probe light from 
the atomic resonance and po(u>) — Vui 2 /2tt 2 is the free 
photon spectral density. The amplitude t%j is a 3 x 3- 
matrix which connects the incoming and outgoing polar- 
izations of the scattering photon. It can be decomposed 
into the sum of a scalar, an antisymmetric and a traceless 
symmetric part. The scalar part is the only one which re- 
mains in the case of the Rayleigh scattering on a classical 
dipole (J = 0, J e = 1). 

In order to characterize the multiple scattering of a 
photon of frequency to, we first calculate its average 
propagator G(u>). The averaging is over the uncorre- 
lated positions r Q of the atoms and over the magnetic 
quantum numbers m Q of the atoms. The first, stan- 
dard average restores the translation invariance. The 



internal average, a trace with a scalar density matrix 
p assuming that the atoms are prepared independently 
and equally in their ground states, restores rotational in- 
variance. The calculation of the average propagator re- 
duces to that of the scalar self-energy = T,(ui)6ij. 
For a dilute enough atomic gas, it can be calculated 
within the self-consistent Born approximation |l]] which 
neglects all possible interference originating from pho- 
ton exchange between distinct atoms: S(w) = NMj t(ui ) 
where Mj = ±(2J e + 1)/(2J + 1) such that M = 1 go|. 
The real and the imaginary parts of the self-energy give 
respectively the photon frequency shift and the elastic 
scattering mean free time r e = — [2Im£(u;)] -1 . The elas- 
tic mean free path l e = r e (remember c = 1) measures the 
average distance between two successive collisions of the 
photons. It is worth emphasizing that l e , up to the scalar 
factor Mj, is independent of the internal quantum struc- 
ture of the atoms. Indeed, the optical theorem assures 
that the total cross section a = l/nl e equals the imagi- 
nary part of the diagonal matrix element of the collision 
amplitude iy(w) which depends on the scalar component 
of this tensor only. 

To go further, we need to calculate the (time inte- 
grated) average probability P(r , r) for a photon to prop- 
agate from a point source ro to a point r. We will calcu- 
late its Fourier transform 
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rewritten as -FW(q) oc G A G R + G A G R Tkk' (q) G A G R . 
The first term is the Drude-Boltzmann contribution. In 
the weak scattering limit kl e 1, the second term 
contains two contributions: the incoherent diffuson 
for which the phase averages to zero, and the coherent 
cooperon T c which still retains an overall phase. In the 
limit of dilute atomic gases, it is possible to write for 
both the diffuson and the cooperon an integral equation 
which generalizes the well-known scalar calculation 
and takes into account the change of the photon polar- 
ization as a result of the existence of internal degrees of 
freedom of the atoms [[ll| . 

In order to describe multiple scattering, we need first to 
determine the atomic intensity vertex U (the differential 
cross section), the average square of the matrix element 
(k'e'|T(w)|ke) = (?' • t(u) ■ e) e l ^ w '> r of the transition 
operator. Averaging over the position and the internal 
quantum numbers yields ||io[| 

U = NTr [p (e 4 .t(uj).e 3 ) (e 2 .t(w).ei)] <5 kl+k3ik2+k4 

= Mj\t(oj)\ 2 eue2je3k£4l Ja,ife£ki+k 3 ,k 2 +k4- (5) 

The vertex U describes the scattering of two incom- 
ing photons (ki,ei) and (k 3 ,e 3 ) into the two outgo- 
ing states (k2,€2) and (k4,C4). The rank-four tensor 
hl,jk — MJ 1 Tr[pdidkdjdi] can be decomposed into its 
three irreducible components Q 
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where the eigenvalues 
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are given in terms of 6j-Wigner symbols. The three basic 
tensors 



T; 



(o) 



1 



(1) 



il,jk 



T; 



(2) 



il,jk 



[Sijdki - SikSji] 



(8) 
(9) 
(10) 



project on the scalar, antisymmetric and symmetric 
traceless components, respectively, of the field intensity 
matrix EiEj. Note that energy conservation imposes 
Ao = 1 for all J, J e . For the pure dipole scatterer 
J = 0, J e = 1, the scattering tensor is proportional to 
the identity, implying that all eigenvalues are Xk = 1. 

Then, the propagation of the intensity between succes- 
sive scattering events is described by the product of the 
average photon propagators, namely 
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where Ay = 5y — r^j describes the transverse projection 
of the polarization vector. In the diffusive limit ql e <C 1, 
the intensity propagator G_gju can also be decomposed 
using the basic tensors d^-^0[) and we recover the eigen- 
values bo — 1, b\ = | and 62 = jq as f° r the classical 
dipolar case 

The summation of the geometric series built from the 
elementary scattering event, i.e. the tensorial product of 
Iil,jk and Gujk, yields in the diffusive limit ql e -c 1 for 
the diffuson 
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where the propagators of the three eigenmodes are 
1 1 
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Here, = /g/3r tr is the diffusion constant, involving the 
resonant transport time scale r tr = l e |l^,|l^]. The 

polarization relaxation times 
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describe the depolarization of the initial light beam and 
are related to the depolarization factors calculated in the 



classical Rayleigh case H|. Among the three modes, 
we have t ( 7 1 (0) = 0, and the corresponding scalar mode 
Ao(g) diverges for small q. It is the singlet Goldstone 
mode associated with the local conservation of the num- 
ber of photons. The antisymmetric mode K = 1 is anal- 
ogous to the triplet mode obtained in the spin-orbit scat- 
tering in electronic systems || . 

A finite time r<j leads to an exponential attenuation 
of incoherent diffusion contribution. In order to identify 
unambiguously a dephasing time, we have to consider the 
cooperon contribution T c . It accounts for the interference 
of amplitudes associated to two time-reversed multiple 
scattering paths of the photons and can be calculated 
along the same lines. The atomic intensity vertex has a 
form similar to (|^), with the new eigenvalues defined in 
terms of 9j-symbols, 
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The cooperon contribution reads 
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where q c = k + k' is the total momentum and the po- 
larization vectors 63 and 64 have been exchanged. The 
propagators of the three eigenmodes are given by 
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The overall depolarization described by the times Td{K) 
affects the cooperon as well. But the cooperon involves 
an additional contribution, the dephasing time 
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This expression constitutes our main result. The limit 
of the classical Rayleigh scattering J = corresponds 
to Xk — Xk — 1 an d thus to the absence of dephasing 
(t>(0) = 00). Of particular interest is the value of 7^(0) 
associated to the intensity of the field. Explicitly, we find 
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An absence of dephasing only occurs for the classical 
dipole J = and in the semi-classical limit J = J e — * 00. 

In experiments on coherent backscattering (CBS) of 
light by atoms without an internal degeneracy (J = 0), 
the optimal CBS enhancement factor is found in the po- 
larization channel of preserved helicity [l5| . The present 
analysis allows to explain the unexpected experimental 
observation that the CBS enhancement factor for cold 
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rubidium atoms ( J = 3, J e = 4) in the channel of pre- 
served helicity is lower than in the channel of flipped he- 
licity ||. Indeed, following O], we define a CBS contrast 

function C{n) — (0)/r|j(0) as the ratio of cooperon 
and diffuson contributions at the nth scattering order. 
For large n, only the largest eigenvalues b = A = 1 
of the diffuson mode expansion will contribute for both 
polarization channels. For the channel of flipped helicity 
(_L), the contrast function behaves as C±(n) ~ 
and decreases exponentially. In the channel of preserved 
helicity (||), all irreducible modes contribute a priori: 
C||(n) ~ Xo — 3(6ixi)" + f (62X2)"- In the case of classical 
dipole scatterers such that xk = 1, the parallel contrast 
stays optimal, Cii(n.) ~ 1. But in the case J = 3, J = 4, 
the product 62X2 = j§ is much larger than Xo = ^ aim 
therefore dominates. The contrast in the perpendicu- 
lar channel then is higher than in the parallel channel, 
C_l/C|| ~ 6, independently of the scattering order n. 

In terms of the above definitions, okXk = 
[1 + T t i/Td(K) + th /T<j,(K)]~ . We can therefore con- 
clude that the loss of CBS contrast in perpendicular po- 
larization channels is not a dephasing process since it 
persists even for J = where t7 (K) = 0, but still 
r d _1 (2) > 0. Since the cooperon and diffuson contribu- 
tions probe different field correlators, it should rather be 
regarded as an effect of polarization decorrelation. On 
the contrary, an atomic internal degeneracy affects the 
contrast in channels of parallel polarizations and there- 
fore appears as a particularly neat realization of a micro- 
scopic dephasing mechanism. 

In summary, we have identified dephasing times in the 
weak localization of light by cold atoms. They depend in 
a remarkably simple manner on the internal atomic de- 
generacy and the field polarization mode. It is of interest 
to notice that the dephasing times can become negative 
(since for instance Xo = — § for J e = J = Although 
this does not change the sign of the contribution to the 
coherent backscattering cone, it might change the sign of 
the correction to the bulk diffusion constant D Jig] , sig- 
nalling an antilocalization contribution like for the spin- 
orbit correction to the conductivity in metals. This may 
lead to an unusual behavior at the Anderson localization 
transition edge. This is quite similar to the case of spin- 
orbit H in disordered metals (symplectic limit) except 
for the fact that here the problem is richer so that the 
critical exponents depend on the Zeeman degeneracy J 
and J e of the atoms |l7j . 
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